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Abstract 

Formulae are provided to express any two-loop scalar integral with 
arbitrary masses and arbitrary external momenta in terms of an inte- 
gral of one fairly simple function or of its derivatives. Such integrals 
can be calculated numerically with high precision. Good agreement is 
found with known analytical expressions of specific two-loop diagrams. 
To prove the effectiveness of these techniques, the two-loop selfenergy 

2 

of the Higgs boson to order {g 2 ^f-) 2 is calculated, and the correctness 

. Tn ' w 

of the result is checked by using the unitarity of the S matrix. This 
result allows one to disentangle the leading corrections to the shape of 
the Higgs resonance, and to derive a perturbative bound of ~ 1.2 TeV 
on the mass of the Higgs boson. 



1 Introduction 



Perturbation theory is practically the only known way to compare field the- 
ory to the experiments. Nonperturbative approaches, like lattice field theory, 
are still far from reaching the precision level of perturbation theory or of the 
LEP experiments. The need for techniques to reliably calculate higher order 
radiative corrections is therefore hard to overemphasize. 

In the case of QCD, two-, three- and even four-loop calculations were 
performed (l|, |||, but the massive case is much more difficult, and an ana- 
lytical solution of the general two-loop integral does not exist. 

If the particles on the external lines are light, so that one can perform 
an expansion around zero in the external momenta, the Feynman integrals 
can be carried out analytically. Such closed formulae made it possible for 
instance to calculate certain two-loop electroweak corrections to the p pa- 
rameter H , to the masses of the vector bosons , and to the selfcouplings 
of the vector bosons ||. 

There are however physically interesting cases when such an expansion 
cannot be performed. One obvious example is the Higgs physics. Knowledge 
of higher order quantum corrections in the Higgs sector is interesting for 
several reasons. First, if the Higgs particle is heavy, they will be numerically 
important and will play a role in the Higgs searches at future colliders. 
Second, as the Higgs mass increases, the Higgs sector becomes strongly 
selfcoupled, and at some point the perturbative expansion will break down. 
Knowledge of the higher order corrections is the only way to find out up 
to which point the calculations by Feynman diagrams still can be trusted. 
Third, having heavy particles on the external legs is the only way to see the 
details of the symmetry breaking mechanism, which are otherwise hidden 
by the screening theorem [||, f?J. 

Considerable effort was devoted lately to the calculation of massive two- 
loop Feynman diagrams. Analytical results exist for certain (mainly selfen- 
ergy) diagrams evaluated at special values of the masses and of the external 
momenta. Small and large momentum expansions were worked out, too 
||] , although the small momentum expansion is a rather trivial problem, 
since the propagators can be expanded in the external momenta, and closed 
expressions for the zero momentum case are available. For the two-loop 
selfenergy diagram with three propagators, an expression in terms of gener- 
alized hypergeometric functions has been found recently [ pX[j , 

These methods refer to special cases and were so far of little use to 
calculate physically relevant quantities. 
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Recently a numerical approach to the so-called selfenergy master dia- 
gram was proposed which reduces this diagram to a two-fold integral [11]. 
It is possible to extend this method to all selfenergy diagrams |12(] . The 
generalization requires a special treatment of different topologies, and also 
numerical problems were reported in the imaginary part for certain values 



of the masses and momenta [12|. 



One would ideally want to use one single method for any topology of 
two-loop diagrams. This would allow one to write computer programs to 
generate and compute automatically all Feynman diagrams needed for a 
certain Green's function. Furthermore, one would like to push the analytical 
integrations as far as possible to increase the accuracy of the final result. 
The Monte-Carlo techniques allow integrations over many variables, but 
their accuracy is usually limited. 

This paper proposes a method to reduce any two-loop massive Feynman 
diagram to a standard form involving essentially one basic function and its 
derivatives, by use of the Feynman parameters. Further, this basic function 
is integrated analytically as far as possible. Its analyticity properties are 
discussed, leading to efficient methods to perform the integration over the 
remaining Feynman parameters numerically. Good agreement is found with 
known analytical results. As an exercise, we calculate in the last section 

2 

the selfenergy of the Higgs boson to order (g 2 —?-) 2 , check the result by 

m w 

using the unitarity of the S matrix, and extract the leading corrections to 
the shape of the Higgs resonance. 



2 Reduction formulae 

In this section we provide formulae to reduce any two-loop scalar diagram 
with arbitrary masses and arbitrary external momenta to a standard form. 
This is an integral over Feynman parameters of one basic function or of 
its derivatives. In one special case (the two-loop integral with three prop- 
agators), a second function needs to be introduced. This second function 
appears rarely in practical calculations and can in some cases be evaluated 
analytically. 

The methods described in this section are easy to encode into an alge- 
braic computer program like FORM [jilj| ], in order to generate and handle 
the large number of Feynman diagrams which may be required to calculate 
Green's functions of physical relevance. 

The most general form of a dimensionally regularized two-loop scalar 
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diagram in Minkowski space is: 



/ = fd n pd n qf\- — - 9 ^ TT 7 ; To 2 

J fj[ W + ki > ~ m i + 111 j=i W + ka +j) - m a+j + %r l 

C 1 

n (r + k a+b+i ) 2 - m 2 a+b+l + i V ' ^ 

where r = p + q is the sum of the loop momenta, and {/ci}i=i,...,a+f>+c are 
external momenta, not all of them independent. 

As an obvious step, it is sometimes possible to reduce the number of 
propagators by partial fractioning: 
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[(p + k) 2 — mf] [(p + k) 2 — m|] 

11 1 
m\ — m\ (p+k) 2 —m 2 (q + k) 2 — m\ 

One then combines all propagators with the same loop momentum by 
using Feynman parameters: 
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[ip + ki) 2 - mf] a i [(p + k 2 ) 2 - m|] Q 2 



r(ai + a 2 ) f 1 ^ x^-^l-x)^- 1 (3) 



T(a l )T(a 2 ) Jo [(p + k) 2 - m 2jai+a 2 

where 

k = xk\ + (1 — x)k 2 
m 2 = xm\ + (1 — x)m 2 — x(l — x)(k\ — k 2 ) 2 ■ (4) 

The integral of eq. 1 can therefore be written in the following form: 

/ = J dX d n pd n qP(xi,x 2 ,- ■ ■ ,x a+ b+ c -3) 

= = ^ = = = , (5) 

[{p + h) 2 - mf\ ai [(q + k 2 ) 2 - ml\ a * [(r + k 3 ) 2 - mj] a 3 
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where P is some polynomial, xi,X2, ■ ■ ■ , x a+ b+ c -3 are the Feynman param- 
eters (we assume that no partial fractioning was possible in eq. 1), and 
X = dx\dx2 ••• cte a +&+c-3- ki, k 2 and k 3 are polynomials of the exter- 
nal momenta and of the Feynman parameters, and m 2 , m?, and 771,3 are 
polynomials of the masses, the Mandelstam variables, and the Feynman 
parameters. 

To simplify the notations, we will drop in the following the tildes, the 
polynomial P, and the integration over the Feynman parameters, keeping in 
mind that the masses and the momenta depend on x±, X2, ■ ■ ■ , x a +6+c-3- 

By redefining the loop momenta, we can rewrite eq. 5 in the following 
form: 



j d n pd n i 



(p 2 - ml) ai (q 2 - m|) Q 2 [( r + fc 3 - fci - k 2 ) 2 - mj}^ 



.(6) 



The Lorentz invariance implies that / depends only on (£3 — k± — A^) 2 . 
It is also obvious that 



J d n pd n q 
J d n pd n q 

I 



[(p + k) 2 - mf] a i (q 2 - m 2 .) 01 * (r 2 - m|) a s 
1 

(p 2 - mj) ai [(q + k) 2 - m 2 ,} * 2 (r 2 - mf) a ^ 



These properties justify the notation: 



G(mi,ai;m 2 ,a2;m 3 ,a 3 ;k ) = 

d >d n q — 2 - — 2 - ^- - — 2 - ^2 ]Q3 , (8) 

where the momenta are now Euclidian. 
We further notice that 



G(mi,ai + l;m 2 ,a 2 ; 7773,03; k 2 ) = 

- — Tr-oG(mi, ai\ 777 2 , ol 2 \ m 3 , a 3 ; k 2 ) 
a\ omf 



(9) 
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Such relations allow one to express all higher order G functions in terms 
of the derivatives of one basic function. We choose this function to be 
G(mi, 2; 1; 1713, 1; k 2 ) , and will denote it in the following by 
Q (mi, 1712,1713; k 2 ). The diagrammatical representation of Q and of its deriva- 
tives is shown in Fig. 1. 

In this way, one case remains uncovered, namely G(m\, 1; 7112, 1; 1713, 1; k 2 ). 
We do not choose this function as the basic function because of its bad ul- 
traviolet properties. Namely, if one naively uses the Feynman parameters 
to calculate it, part of the ultraviolet divergencies will be transferred from 
the radial integration to the integration over the Feynman parameters. 

To calculate G(m\, 1; mi-, 1; ms, l;k 2 ), we start by using the "partial p" 
operation [O]: 



G(mi, l;m 2 , l;m3, 1; k ) = 

d n pd n q ' 



(p 2 + m\)[(q + k) 2 + m 2 ](r 2 + TO3) 

,2 



2 f d n pd n q{ '' 



n J (p 2 + m^) 2 [(g + &) 2 + w- 2 ] (r 2 + m 2 ) 

+ P 2 + p - g | (10) 

(p 2 + m\)[(q + /c) 2 + m^](r 2 + m§) 2 

After some algebra and explicit symmetrization upon m\ *-* 1712 , one 
ends up with the following reduction formula: 



G(m 1 ,l;m 2 ,l;m 3 ,l;k 2 ) = 

-^—{m 2 g( mi ,m 2 ,m 3 ; k 2 ) + m 2 g(m 2 , mi ,m 3 ; k 2 ) 
6 — n 

+mlg(m 3 ,mi,m2]k 2 ) + T(mi,m2,mz]k 2 )} , (11) 

where 



f(mi,m 2 ,m3;fc 2 ) = 

d n pd n ( P + q).k 

(p 2 + 7nf)[(q + k) 2 + ml](r 2 + m|) 2 

The function .F must be symmetrical under permuting m\, 1712 and 
1713 due to the original symmetry of G(mi, 1; m2, 1; 7713, 1; A; ) and to the 
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77i2 7713 symmetry of £/(mi, 7712, 7713; k 2 ). Another property of T is that it 
vanishes in the limit k 2 — ► 0. Indeed, the function (mi, 7712,7713; A; 2 ) reduces 
in this limit to the symbol (mi, mi |m2| 777,3) of ref. ||. According to ref. 
the following formula holds: 



G(mi,l;m 2 ,l;m 3 ,l;0) = 

{m\ g(mi,m 2 ,m 3 ; 0) + m\ g(m 2 ,m 1 ,m 3 ; 0) 

3 — n 

+m|^(m 3 ,mi,m 2 ;0)} . (13) 

This proves the statement about the limit of T . 

We can at this point express any two-loop Feynman diagram as an inte- 
gral over Q functions or over the derivatives of Q. In the case of the two-loop 
diagram with three propagators and nonvanishing external momentum, one 
T function appears, too. 



3 Analytical integration 

In this section we integrate the T and Q functions analytically as far as 
possible. One will be left in the general case with an integral over one 
Feynman parameter. Only in special cases can this last integral be solved 
analytically. However, a careful examination of the integrand's analyticity 
properties allows one to write fast computer programs to evaluate these 
functions numerically with high precision. 

We start with the Q function. After introducing two Feynman parame- 
ters and integrating over the loop variables, one obtains: 



g(m 1 ,m 2 ,m 3 ; k 2 



f d n pd n q - 

J (p 2 + m\ ) 2 \(q + k) 2 + ml] [(p + q) 2 + m|] 



r(3-§)J 



TV 9 - Vl\ rl 
^(vrm 2 )™- 4 M 2 > 



{r(5-n). 



dx 



dy[x{l-x)]2 2 y{l-yf 



y 2 n 2 + fi 2 



[y(l-y)K 2 + y + (l-y) f i 2 ]5-n 



[y(l - y)n 2 +y + (l-y)[i 



(14) 
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where the following notations were introduced: 



9 ax + 6(1 — x) 
x(l — X) 

m 2 . _ m l 



mf mf mf 

We then expand in e = re — 4 and integrate over y . After some tedious 
but elementary algebra, one obtains: 



G(m 1 ,m 2 ,m 3 ;k 2 ) = 

^{ 4 + -I" 1 + 2 7 + 21og(7rm?)] + ~ + ^ 
e z e 4 12 

+ i[-l + 27 + 21og(7rm?)] 2 -l+ 5 (m 1 ,m 2 ,m 3 ;A; 2 )} , (16) 



where 



1,1 1 1 
g(mx,m2,m 3 ;k z ) = I dx[Sp(- ) + Sp(- 



o 1-2/2 

+2/1 log ^ : + 2/2 log V2 : ] , (17) 
2/1-1 2/2-1 



and 



1 + K 2 - /i 2 ± y/A 



Vl > 2 ~ 2^ 



A = (1 + k — A* ) 4k — Aik i] . (18) 

In general the finite part 5 cannot be integrated analytically. However, 
this is possible in the trivial case when k 2 — ► . One can convince oneself 
that eq. 17 reduces in this limit to 

f 1 /u 2 

g(mi,m 2 ,m 3 ;0) = dx [ 1 + Sp(l - fJ 2 ) - - nlogAi 2 ] , (19) 

Jo 1 — jtr 

which is essentially the finite part of the function (mj, mi\m2\m 3 ) of ref. ||. 
In this case the x integration can be carried out, and gives: 
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g(m 1 ,m 2 ,m 3 ;0) = 

1 - ^ log a log 6- a + r— 1 [Sp(-—) + Sp(-—) 
2 vA/ v\ ui 



1 r, Uo 1 , 9 V2 1 , ,1*1 1 , n 1i 2 7T 2 . , . 

+ log 2 ^ + -log 2 ^ + -log 2 ^--log 2 ^ + — , 20 
4 fi 4 ui 4 «i 4 



where 



ui,2 = + 

A' = l-2(a + b) + (a-b) 2 . (21) 

This agrees with the results of ref. || . The result of eq. 20 implies that 
any massive two-loop Feynman diagram can be calculated analytically if all 
Mandelstam variables vanish. 

Let us now turn to the T function. Eq. 12 can be rewritten in the 
following form: 

T{m\,m2,m?,\k 1 ) = k 2 G(mi,m 2 ,m 3 ;k 2 ) + J r '(mi,m 2 ,m 3 ;k 2 ) 

F'( mi ,m 2 ,m 3 ; k 2 ) = 

d n d n P-fc (22) 

P q [(p + k) 2 + ml} 2 (q 2 + ml) [(p + q) 2 + m 2 ,] 

We combine the propagators by using Feynman parameters, and inte- 
grate over the loop momenta: 



J r '{mi,m 2 ,m 3 ; k z 



-fcV(™ 2 r 4 ^y| £ dx £ dy [x(l - xp- 2 y 2 (l - y) 2 -f 
{ r(5-n) y2K2 + » 2 



[y(l - y)K? + y + (1 - y)^ 2 f- n 
n„,. x 1+ 2 



2 [y(l - yjK 2 + y + (1 - y)n 2 f n 
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with the notations of eq. 15. 

One next expands in e = n — 4 and integrates y out. Using the result for 
6(mi,m2,m3;fe 2 ) , one finds: 



f(mi,m2,m3;fc 2 ) = 



A ; 2 7r 4 {-l + ^-^[7 + log(7rm?)]+/(m 1 ,m 2 ,m 3 ;A ; 2 )} , (24) 



where 



/(mi,m 2 ,m 3 ; k 2 ) 




yi - 1 




2/2-1 



(25) 



As already mentioned, the T function must be invariant under permu- 
tations of the masses, and has to vanish in the k 2 — > limit. While the 
m2 <-> ma symmetry is obvious, the mi <-> m2 symmetry is less trivial, and 
was checked numerically. 

The / function is needed only for the two-loop diagrams with three 
propagators, and can be integrated analytically for special choices of the 
arguments. 

4 Analyticity and numerical integration 

We found that the finite part of the basic function Q is expressed in the 
general case through a one dimensional integral: 




g(m 1 ,m 2 ,m 3 ;k 2 ;x) 



Sp( 



1 



) + s P ( 



1 



) + 



1 -y\ 



1-2/2 



2/1 log 




+ 2/2 log 



2/2-1 



1 + K 2 - n 2 ±\fA 

L ' 2 - 2^ 

A = (1 + k 2 - fi 2 ) 2 + 4kV - AiK 2 7] . 



(26) 
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It is straightforward to calculate the derivatives of g with respect to 777 1 , 
7712 and 7773, which are needed to evaluate Feynman diagrams with more 
than four propagators. 

In eq. 26, k 2 > corresponds to Euclidian momenta. When one rotates 
k back to Minkowski space to calculate the physical Green's functions, the 
function g and its derivatives develop an imaginary part above the threshold 
—k 2 > (mi + m2 + 777.3) 2 • This is readily seen from the diagrammatical 
representation of fig. 1. 

The typical behaviour of g as a function of the integration variable x 
is shown in fig. 2. It has singularities of the logarithmic type at x = and 
x = 1 . Once they are extracted through a suitable change of variable, the 
function g can easily be integrated numerically. 

Problems start to appear when trying to evaluate Feynman diagrams 
with more than 5 propagators, above the threshold. For the case of six 
propagators there are two Dirac 5 functions appearing at the points where 
the derivative of the function g is discontinuous, and starting with seven 
propagators one deals already with the derivatives of the 5 distribution. 

To circumvent this kind of problems, one has to analytically continue 
the g function in the complex plane of the Feynman parameter x, and to 
choose an integration path which avoids the singularities. 

The singular points of g, other than x = and x = 1, are given by the 
roots of the discriminant in eq. 26. The equation A = has four roots: 

xi,2 = ^2 [" a + 6 + tJ>i±>J(a,-b- tf) 2 -46//?] 

x 3A = 2j?[-a + b + f£±y/(a-b- n\) 2 - 46/x 2 . ] , 

[i\ 2 = 1-k 2 t2V^ • (27) 

These are branching points for each separate term in the expression 26 
of g . However, the g function itself has only two branching points at x\ 
and X2, because the singularities at X3 and X4 are compensating among the 
four terms of g in eq. 26. 

This behaviour is related to the causality of the Green's functions. The 
causality condition is expressed by the ir\ prescription in the Feynman prop- 
agator, wich means to replace all masses m 2 by m 2 — irj. An equivalent way 
to impose causality is to calculate the Euclidian Green's functions and to go 
afterwards to physical momenta, approaching the cut on the positive real 
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axis from above. This amounts to making the replacement k 2 — > —k 2 — irj in 
eq. 26. These two prescriptions ought to be equivalent, and therefore have 
to fix the location of the physical singularities with respect to the real axis 
in the same way. For x\ and x%, at — k 2 > 1 both prescriptions lead to the 
same change, and therefore these are the singularities of the g function. For 
X3 and X4, the two prescriptions would lead to opposite changes in the imag- 
inary direction. Since causality fixes the location of the singularities of the 
Green's function uniquely, X3 and X4 cannot correspond to real singularities 
of g. Therefore g is analytical at these two points. 

Fig. 3 shows the typical behaviour of g above the threshold as a function 
of the complex Feynman parameter x . One can now choose a complex 
integration path on the physical Riemann sheet along which the integrand 
is free of singularities. On such a path the g function is analytical, and so 
are its derivatives of any order, which are requested to calculate diagrams 
with more than four propagators. 

To calculate Feynman diagrams, one needs to integrate the g function 
or the derivatives of g over the Feynman parameters which were introduced 
to combine the propagators with the same loop momentum. In general the 
arguments of g will take negative values, too, corresponding to imaginary 
masses in eq. 26. This corresponds to various thresholds of the diagram. 
The points where the mass arguments vanish must also be avoided, since g 
and its derivatives may have mass singularities at these points. One needs to 
make sure the integrand remains always on the physical sheet. Therefore one 
always chooses the integration path according to the Feynman prescription 
for the masses of the propagators: m 2 — ► m 2 — irj . 

A computer program based on this method was written and checked 
upon known analytical and numerical results for specific diagrams. 

The topology of fig. 5 a) in the case s = 1, m 2 = m 2 = . . . = m 2 — > , 
which is related to the large momentum limit, was calculated and found to 
agree with the known value of 6£(3) |14| |. Also the case s = 1, m 2 = m 2 = 



771,4 = 1) m 2 = m 5 — * agrees with the analytical result ir log 2 — §C(3) pd[| . 
For the diagram of fig. 5 h) , an analytical formula in terms of Lauricella 



functions was recently derived |i0| . The authors of ref. [10] calculate an 
ultraviolet finite combination Ti23at of four such diagrams. T123N can be 
expressed in terms of 12 g functions and 4 / functions. Perfect agreement 
was found with the numerical values which are given in ref. for a range 
of masses and momenta. 

Finally, some comments on the numerical integration are in order. The 
g function can easily be integrated through virtually any algorithm after 
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extracting the singularities at the ends of the integration path. These are 
of logarithmic type, and can be extracted minimally through a change of 
variables of the type: 

t = x(24 - 24 log x + 12 log 2 x-A log 3 x + log 4 x) . 

A simpler nonminimal change of variable, like t = \fx , can be more 
handy. It puts the integrand to zero at the ends of the integration path. 

Because the integrand is smooth along the integration path and free of 
violent variations, the g function and its derivatives can be integrated nu- 
merically very fast, using a small number of points. Typically some 120 
points were requested to calculate these functions to more than 8 digits. 
This takes about 50 ms on an IBM RISC 6000 workstation. Integrating the 
g function or its derivatives further in order to calculate Feynman diagrams 
with more propagators poses also no problem as long as the appropriate 
complex integration path is used in order to avoid the singularities. This 
yields also a smooth function which can be integrated with high precision. 
For less than four Feynman parameters, an adaptative deterministic algo- 
rithm was used, yielding accurate results already with a small number of 
points. For more than four dimensions the Monte-Carlo techniques become 
superior, but one cannot hope to obtain an accuracy better than 3-4 digits. 



5 The selfenergy of the Higgs boson 

To show how the techniques of the previous sections work, and to prove that 
they can be used to perform reliable and accurate calculations of physical 
relevance, we calculate in this section the selfenergy of the Higgs boson at 

2 

order (g 2 —^) 2 - We first calculate the on-shell selfenergy, check the cor- 

m w 

rectness of the results by using the unitarity of the S matrix, and extract 
the mass counterterm. We then calculate the momentum dependence of the 
Higgs selfenergy to see its analytic structure, and check some of its asymp- 
totic properties. The result allows one to extract the leading corrections to 
the shape of the Higgs resonance. This is an effect which does not appear in 
the selfenergy of the Higgs boson at one-loop level. The leading corrections 
to the Higgs shape become large if the Higgs mass is of the order of 1.2 TeV, 
indicating a strongly interacting theory, in agreement with well-known Born 



level and one-loop results fig, |19[ . 

Since we are interested in the leading contribution in rriH, the most 
natural choice is to work in Landau gauge. In this gauge only the Higgs 
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sector survives at leading order in nin- The diagrams containing gauge, 

2 

fermion, or Fadeev-Popov fields do not give contributions of order {g 2 —?r) 2 

m w 

in this gauge. 

However, to avoid problems with massless Feynman diagrams which can 
be traced back to the arbitrariness of / d n p ^ within the dimensional regu- 
larization, we choose to keep a small gauge parameter £ during the compu- 
tation, and take the limit £ — > in the final result. This amounts to keeping 
a finite mass of the Goldstone bosons. The diagrams involving gauge and 
Fadeev-Popov fields can safely be omitted, since they do not generate finite 
terms in the limit £ — ► 0. This provides also a useful check of the compu- 
tation, since all poles and logarithms of the gauge parameter have to cancel 
in the final result. 

The Lagrangian of the Higgs-Goldstone system is: 

C = ^H )(d»H ) + ±(d tl z )(d»z ) + (d IM w+)(d»w r) 

-9 —2-^[w^w +-z + -H + ——H + ] (28) 

with obvious notations. St is the tadpole counterterm needed to ensure that 
< 0|iT|0 >= v at all orders. 

The renormalization is performed by splitting the bare quantities in 
renormalized quantities and counterterms: 



H = z]j 2 H 
zo = Zq z 

wo = Zq w 

2 2 x 2 

m Ho = m H -Sm H 

m w = m w ~ ^ m w ■ ( 29 ) 

Note that the field renormalization constants of the charged and the 
neutral Goldstone bosons can be chosen equal due to the remnant 0(3) 
symmetry of the Lagrangian. 



We adopt the tadpole renormalization strategy described by Taylor [16]. 



The condition that the vacuum expectation value of the Higgs field should 
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not receive quantum corrections amounts to neglecting the tadpole insertions 
altogether and subtracting the value of the Goldstone bosons selfenergy at 
zero momentum from all scalar selfenergies. 

The counterterms are fixed uniquely through the renormalization condi- 
tions. We use an on-shell renormalization scheme with field renormalization, 
and the physical masses of the Higgs and of the W bosons as input param- 
eters. Note that the gauge coupling constant does not get renormalized at 

leading order in g 2 ^-. 

m w 

By evaluating the one-loop selfenergies of the Higgs, Goldstone and W 

2 

bosons, one finds the following renormalization constants at 0(g 2 ^^): 



xdi-ioop) _ 2^|l i m H ™?h s _ A , 3 _ 3_7 
° ~ 9 ~ 2 K *—*> 16^ 1 4e + 8 8 




5mT l ° OP) = ,4(4) £/2 ^{--3+^ + ^ 
H m w 4^ 16vr 2 e 2 8 

ro 3 7 , 3 7 2 3^/3^ 3\/~3~/ir 

+e[3 -^ + ~8 ~ + ^^~ 

tt 2 3vg_C7(|) 3V37T log(3) 
~16~ 4 + 16 " 

, 2(l-Ioop) _ 2 ™ff / "4 W2 4rl,/ 3 , 7 \ , 

dm w ~ 9—^4^2) 16^ [ 8 + U 32 + 16/ J 



w 



c 7 (l-loop) J™>H (™H_\e/ 2 1 r 3 TTV3 f 3 3 7 

6Zh - 9 ^ { ^> T6^ ! 2"^ + e[ "2 + T 

3V3tt v^tt vg_C7(|) V37TlQg(3) 
+_ 16 8 + 2 8 " 

6Zg ~ l6^ i_ 8 +e V 32"l6j j • (30) 

One needs the one-loop counterterms at 0{e) because they combine with 
the ^ terms at two-loop order to give finite contributions. 

We can now turn to the actual two-loop calculation. 

The two-loop tadpole counterterm can be calculated from the diagrams 
of fig. 4. Alternatively, one can calculate the w oi z selfenergies at zero 
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momentum. In any case, the tadpole counterterm can be evaluated analyti- 
cally, since it involves only the g(rni, fn%, 7743; 0) function and its derivatives. 
It was checked that all three ways to calculate the tadpole counterterm lead 
to the same result: 



c.(2-loop) _ I 2 m H \2 1 m H \e m H r ^5 , \ i ^3 

{9 m 2 / l 4V j (16vr 2 ) 2 [ 16e 2 + e [ 8 

457 9y/3ir, 609 337 457 2 45V3~tt 

H H ) H H 

16 16 ; 128 8 32 64 

9V3 7 7r 3^ 2 2U/3C7(|) 9V37rlog(3) 

16 32 32 32 U ' 

Note that this expression disagrees with the results of ref. [|l7|| p|. 
The main task is to calculate the two-loop diagrams of the Higgs selfen- 
ergy. The topologies involved are shown in fig. 5. 

Subtracting the tadpole counterterm from the diagrams of fig. 5, and 

performing the necessary integrations numerically, one obtains the following 

2 

Higgs selfenergy at 0((g 2 ^-) 2 ): 



2 rn^ 2 m|_ rr^_ _9 + 3 _ _ ^ 

(1d7t z J z e z die 
-4.6298(7) -i 0.4124(5)] . (32) 



The real part of this expression gives the two-loop Higgs mass countert- 
q. The absor 
decay width JO 



2 

erm. The absorptive part agrees with the 0{g 2 —^) corrections to the Higgs 



T(H -» W + W~,ZZ) = 

» 2 4 T |-[ 1+ '' 2 4A(^-^ + ^)] • (33) 
m^r 128 tt my V 8ir z 16 8 48 

1 At least two mistakes exist in ref. j^]. Their diagram £20 has the value and 

not C(2) + 4 log 2; and £3 = — ^- + 2y / 3C ? /(f ), so it does not contain the logarithm of 
the golden ratio, as claimed in those papers. These two diagrams enter the expressions of 
both Goldstone and Higgs selfenergies, and hence affect all their results. 
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Note that this is not only a check of the imaginary part of the diagrams. 
The path in the Feynman parameters integrals was varied in large limits 
to check that the result is the same. This proves that the integrand is 
indeed analytical. Since an analytical function is uniquely determined by 
its imaginary part, it is to be concluded that the real part of the numerical 
integration is correct as well. 

The momentum dependent selfenergy of the Higgs field is given by: 

X%- loop \k 2 ) = 



(9 2 ^f [-4 + ^"(168 



m l w ' Mtt^ 2/ (16vr 2 ) 2 L e 2 32e m 2 H 

k 2 

-967-24^) + £ / (_)] . (34) 



The behaviour of its finite part S^(-t-) is shown in fig. 6. It has the ex- 
pected analytical structure. £/, as well as the one-loop selfenergy, displays 
unphysical singularities at k 2 = and k 2 = 4m 2 H , which cancel in the full 
Green's functions. The imaginary part of the selfenergy grows like k 2 at 
large momenta. Considering the asymptotic behaviour of the phase space 
factor of the three-body decay, this agrees with the Cutosky rule. 

The real part of the derivative of the selfenergy at k 2 = m 2 H is absorbed 
in the wave function renormalization of the Higgs field. On the contrary, its 
imaginary part is physical, and gives the leading radiative corrections to the 
Higgs shape. The momentum dependence of Im T,f(k 2 ) describes corrections 
to the Breit-Wigner shape, corresponding to an energy dependent width. 
By expanding the selfenergy around the peak and keeping only the first 
derivative, one obtains the following correction to the Higgs propagator: 

1 1 



k 2 — m 2 H + imnTn k 2 — m 2 H + imuTu + i(k 2 — m H)^'H ' 



16 7r 2 m 2 ^ dk 2 



1-0 J 2 2 < 35 

16 7r 2 m^y 



One readily sees on dimensional grounds that such an effect does not 
appear at the one-loop level, because the phase space factor is s independent, 

2 

and therefore no s-dependent width can occur at order g 2 —*r . 
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mH 
[UeVJ 


peak 
snirt 
[UeVJ 


peak 
snitt 

[/o\ 


height 
increase 

L/oj 


400 


-.072 


-.018 


.0020 


500 


-.35 


-.070 


.012 


600 


-1.3 


-.21 


.051 


1 uu 


'i s 
-o.o 


-.04 


1 8 
. 15 


800 


-9.9 


-1.2 


.51 


900 


-23 


-2.6 


1.3 


1000 


-49 


-4.9 


3.0 


1100 


-97 


-8.8 


6.5 


1200 


-178 


-15 


13 


1300 


-305 


-23 


25 


1400 


-487 


-35 


45 


1500 


-721 


-48 


78 


1600 


-997 


-62 


130 



Table 1: Leading corrections to the shape of the Higgs resonance 
1-5 -i 1 77 „ | 2 . The Breit-Wigner resonance with constant width calcu- 

1 p — m H — ilmlj 1 ° 

2 

lated at order (g 2 —/r) 2 is compared to the resonance corrected for the energy 

m w 

dependence of the two-loop selfenergy. 



The numerical value of T' H given in eq. 35, obtained from the momentum 
dependence of the Higgs selfenergy shown in fig.6, agrees with the exact 
result which can be derived by using the Cutkosky rule: 



- 2 

(.9 



2 m H \2 




1.002245142 —iLy -— — . (36) 



This correction shifts the peak of the Higgs resonance towards a lower 
energy, and increases its height a little, as shown in fig.7. Further corrections, 
which come from the second derivative of the Higgs selfenergy, are much 
smaller, and do not affect the position of the resonance — they only make 
it deviate more from the Breit-Wigner shape. 
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The magnitude of these corrections for various values of m# is given in 
table 1. The corrections become large if the Higgs boson is heavier than 
~1.2 TeV, signaling strong couplings in the electroweak symmetry breaking 
sector. This agrees with well-known results based on the unitarity violation 
in vector boson scattering at tree level [18], and on the magnitude of one- 
loop radiative corrections to the Higgs width [19|. On the other hand, this 
perturbative bound is considerably lower than the value of 3 — 4 TeV derived 
from two-loop corrections to the p parameter and to the selfcouplings of the 
vector bosons [|,§. This was to be expected, since the latter are corrections 
to low energy parameters, therefore subject to the screening theorem — the 
leading contributions in run cancel. 



6 Conclusions 

We described a method to calculate two-loop massive Feynman diagrams 
which can be applied, at least in principle, to any diagram. 

The diagrams are integrated analytically as far as possible. The neces- 
sary formulae are easy to encode into an algebraic computer program. The 
remaining integrals are made smooth and free of singularities by appropriate 
changes of variables and convenient choices of the integration paths. The 
resulting formulae are suitable for numerical integration since they have a 
small variance. 

The techniques were compared to analytical and numerical results for 
certain two-loop diagrams. In all cases good agreement with the known 
results was found. 

The complexity of Green's functions which can be calculated by using 
this method is limited only by the time needed to perform the numerical 
integrations. 

For less than four Feynman parameters, deterministic numerical inte- 
gration methods were used, which lead to fast and accurate answers. In 
particular, two and three point functions appear to pose no problem. 

For more than four Feynman parameters to integrate over, the Monte- 
Carlo integration is expected to be faster, but one has to content oneself 
with 3-4 digits accuracy. The time needed to calculate the integral is not 
expected to grow too fast with the number of Feynman parameters, since the 
speed of the Monte-Carlo integrations depends primarily on the variance of 
the integrand, and not on the dimension of the integral. 

To show that the method can be used to calculate Green's functions 
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of physical relevance, which usually involve many diagrams and large mass 
splittings, it was used to calculate the selfenergy of the Higgs boson and to 
extract the leading corrections to the shape of the Higgs resonance. 

The corrections to the Higgs shape imply a shift of the resonance towards 
a lower energy. In the range up to 800 GeV, which will be covered by the 
four lepton events at LHC, the shift is quite marginal, at best at 1.2% level. 
The shift becomes increasingly important for heavier Higgs bosons, which 
can presumably be searched for at LHC by looking at the jet decay modes. 
However, to determine the implications of the peak shift for the search for 
a very heavy Higgs boson would request a full analysis of the production 
mechanism. 

The shape corrections grow with ran, becoming rather large for mu ~1.2 
TeV, where the perturbative approach eventually breaks down. 
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Figure captions 



Fig.l The diagrammatical representation of the function £?(mi, r?i2, 777,3; k 2 ) 
(a) and of its higher order derivatives (b). 

Fig. 2 The behaviour of g(mi,m2,m3; k 2 5 OC^ ctS Si function of x under the 
threshold (a) and above (b). The solid line represents the real part, and the 
dashed one — the imaginary part. 

Fig. 3 The real (a) and the imaginary (b) parts of g(mi,m2,ms; k 2 ;x) as a 
function of the complex Feynman parameter x. The logarithmic singularities 
at x = and x = 1, as well as the two branching points close to the real 
axis between and 1 can be seen. The integration path with the ends at 

and 1 avoids the first singularity by going through the positive imaginary 
half plane, then changes to the negative imaginary half plane to avoid the 
second singularity. 

Fig. 4 The topologies of the two-loop tadpole diagrams. 

Fig. 5 The topologies of the two-loop selfenergy diagrams. 

Fig. 6 The real (solid line) and the imaginary (dashed line) parts of the 
finite part £/• (-t-) of the selfenergy of the Higgs boson. 

Fig. 7 The leading corrections to the shape of the Higgs resonance 
1-2 i 1 .j ^ \ 2 . The Breit-Wigner resonance with constant width (solid 

1 p^— m H — ilmX, 1 ° v 

line), and the correction for the energy dependence of the two-loop self- 
energy (dashed line). 
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